A comment on "Amplification of endpoint structure for new particle 

mass measurement at the LHC" 
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Abstract 

We present a comment on the kinematic variable mcT2 recently proposed in [l|. The variable is designed 
to be applied to models such as i?-parity conserving Supersymmetry (SUSY) when there is pair production of 
new heavy particles each of which decays to a single massless visible and a massive invisible component. It was 
proposed in [l|] that a measurement of the peak of the mcT2 distribution could be used to precisely constrain 
the masses of the SUSY particles. We show that when Standard Model backgrounds are included in simulations, 
the sensitivity of the mcT2 variable to the SUSY particle masses is more seriously impacted for mcT2 than for 
other previously proposed variables. 



If new physics is discovered at the Large Hadron Col- 
lider (LHC) |2j we will wish to measure the masses of any 
new particles discovered. Several methods of measuring 
the masses of such new particles have been suggested (for 
a recent review see 3 ) . 

Here we comment on the mass measurement variable 
mcT2 recently proposed by Cho et al. The variable 
introduced in that paper was employed to constrain the 
participating particles' masses for events characterised 
by pair production of identical-mass heavy particles, Y, 
each of which decays to a massive invisible Xidf) an d a 
single massless visible v(p) daughter, where the symbols 
in brackets label the momenta. 

The experimental signature for this process 



This variable is a novel amalgam of two methods previ- 
ously described in the literature. The procedure of min- 
imising the larger of two quantities — over all partitions 
of the invisible particles' transverse momenta consistent 
with the missing transverse momentum p T — has its ori- 
gins in the 'Stransverse mass' of 0,0] ■ However, whereas 
[J, |i| take the larger of the two transverse masses mr for 
each Y decay, in ([2]) the quantities being evaluated at 
each possible partition of f T are the contransverse mass 
functions mcT- These are defined (for m v = 0) bjQ 
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(pr,t T ,x) = X 2 + 2E(p T ,0)E($ T ,x) + 2p T -$ T , (3) 
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where the transverse energy of a particle with trans- 
verse momentum pt and mass m is given by E(px, m) = 
\J p\ + to 2 , and x is a trial value for the (a priori un- 
known) mass of the invisible particle. The + sign in front 
of the inner product of the momenta in ([3]) distinguishes 
mcT from the usual transverse mass. 

The already bloated dictionary of transverse mass 
variable names is stretched nearer to (or perhaps be- 
yond) breaking point by the addition of the term 
constransverse mass for the quantity defined in ([2]). 

It was suggested in [l[ that by measuring the kine- 
matic end-point of the mcTi distribution it should be 
possible to determine the mass of the parent particle Y 
accurately. This assertion is based on the observation [l| 

1 mcT was originally proposed in Q for the case where both momenta correspond to visible particles and was motivated by particular 
invariancc properties under back-to-back boosts in the transverse plane. The use of mcT for the case where one input corresponds to a 
visible particle, but the other represents the hypothesised momentum of an invisible particle was an innovation of [ij . 

2 Colloquially known as the condition of no 'upstream transverse momentum'. 



is therefore two visible particles (in this paper we assume 
these to be jets and those jets to be massless) and large 
missing momentum. 

Cho et al. suggest constraining the masses by mea- 
suring the distribution of the variable [l| 
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that under the conditional 
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the distribution of tticti is highly peaked at its maximal 
value (under that same condition), 
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2(E T (p ,x)\pa\ - \po\ 
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where \po\ is the absolute momentum of the daughter 
particle in the rest frame of the parent^. 

The central observation of [JJ is that the vhcti distri- 
bution has a sharply defined Jacobian peak at its kine- 
matic cndpoint (|5|) , so a good measurement of that end- 
point position could provide a precise constraint on the 
parent particle mass. This statement is founded on the 
reasonable expectation that the systematic uncertainty 
in fitting a sharp peak ought to be smaller than in fitting 
other distributions which tend to have smaller numbers 
of events near their kinematic endpoints. 

We perform simulations similar to those in [l[ but 
including the most important Standard Model back- 
grounds. We find that those backgrounds are peaked 
in the same region as the signal, so play a much more 
significant role for the process of interest (HJ than might 
be inferred from [l[ . The end-point value (O is relatively 
insensitive to the physical parameter \p 1 , so a rather pre- 
cise measurement of the peak position would be required 
to constrain the parent particle masses. We suggest that 
fitting the endpoint position with the required precision 
is likely to be difficult when the systematic effect from 
uncertain residual Standard Model backgrounds is taken 
into account. 

As is highlighted in Q, a judicious choice of \ m 
needed if the distribution of mcTi is to be sensitive to 
the particle masses. Cho et al. show that if one selects 
a value % ^ \Po\ one loses the sharp peak in the mcT2 
distribution. 

The variation of the endpoint position ([5]) with re- 
spect to \po\ is 
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0|po| 4m%™E T (p , X )\po\ 2 ' 
so if one chooses the other extreme with \ <C |po 1 > then 
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which has very limited sensitivity to the physical param- 
eter \p 1 . Therefore the value of \ that ought to be cho- 
sen should be close to \p \ (rather than for example the 



true invisible particle mass). For x/\Po\ = {0-5, f, 1.5}, 
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■ « {0.02, 0.1, 0.2} respectively. As was recognised 
even with a well-tuned % the peak position does 
not vary much with \po | so an accurate and precise mea- 
surement of tUq^ is required to constrain \po | (and from 
it the particle masses). 

To study the effect that Standard Model backgrounds 
might have on the method proposed, we use a simulation 
similar to that described in Q including a parameterised 
detector response typical of a general-purpose LHC de- 
tector. Herwig++ 2.4.2 is used to produce samples 
of the following Standard Model backgrounds; QCD, ti, 
W — > Z^+jets, Z —> Z + Z~+jets and Z —> w+jets for 
proton-proton collisions at centre of mass energy 7 TcV. 
For a signal we use inclusive production of the SPSla 
supersymmetry benchmark point [Io| with the spectrum 



and decay table calculated by SPheno 2.2.3 [ll( . Jets 
are formed and smeared following the same procedure as 
described in 0} 

The expected peaking of the mcT2(x = 300 GeV) 



variable can seen in Fig. 1(a) after applying only mini- 
mal selection cuts requiring at least two jets with px > 
50 GeV and with pseudorapidity, |t;| < 2.5. We have 
chosen x = 300 GeV, close to po (which ranges from 
240 GeV < \po | < 300 GeV depending on which squark 
mass is used). In the region of the SUSY signal, it can 
be seen that the Standard Model backgrounds are large. 

It can also be seen from Fig. 1(a) that the dependence 
of the position of the peak in the mcT2 distribution on 
the physical parameter \pq | that one is trying to measure 
is very small. This is shown by the different coloured 
arrows on that plot indicating expected signal peak po- 
sitions corresponding to different parent particle masses. 
Even with large variations in the parent particle mass 
(up to 100%) little variation in the peak position is seen. 
This means that a very precise determination of the peak 
position would be required in order to get sensitivity to 
the parent particle mass0 

In Fig. 1(b) we show the mcT2(x = 300 GeV) dis- 



tribution after the application of more selective SUSY 
cuts [TJI proposed by ATLAS. The backgrounds appear 
reduced but are still significant in the signal region, and 
they have a peaked structure similar to the signal. These 
backgrounds will therefore have to be well-understood for 
the successful extraction of the signal shape. 

As was pointed out in [TJ it is possible to sharpen the 
end-point structure of the mcT2 distribution by selecting 
a subset of events for which the 'upstream momentum' 
is small; i.e. those for which (j4]) is approximately satis- 
fied. In Fig. 1(c) we show the mcT2 distribution after 



3 For a two body decay |po| 
SUSY particles respectively. 



2m ™ where my uo and m^ ne are the true masses of the parent and daughter 



This is also true for other values of x- 
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an additional selection requiring \S\ < 30 GcV where the 

upstream momentum S = —f T — Pj< — Pp ■ The incT2 
cndpoint position becomes better defined when this ad- 
ditional requirement is applied, but at the cost of a factor 
of ~ 10 in number of events, and even then the residual 
backgrounds are still not negligible (given the precision 
in the endpoint position that is required). 

To investigate the effect of the Standard Model back- 
ground on the precision with which |po| can be deter- 
mined we parameterised separately the signal s and the 



background b distributions from Fig. 1(c) by two-part 
piecewise Gaussian functions, 



G(^\af):(x> M s < b ); G(//' b , 



5 ) : (x < M s ' b ). 



Even with a high cross section SUSY model (such as 
the SPSla point shown) and assuming perfect knowl- 
edge of the shape of the background^, a significant in- 
crease in the statistical uncertainty is found when the 
backgrounds are introduced. We also investigated the 
case when the SUSY cross-section is a fraction of that 
shown. For lib" 1 the statistical precision on |po| was 
{(±6.8,±9.6), (±8.0,±12), (±14, ±22)} GeV for the cases 
of: original cross-section without and with backgrounds 
included, half of the original cross-section without and 
with backgrounds included, and a quarter of the original 
cross-section without and with backgrounds included re- 
spectively. Uncertainties in the shape of the background 
contribution will further increase the detrimental effect 
of the backgrounds. 

These difficulties with the constransverse mass vari- 
able are not shared by other variables which have pre- 
viously been proposed for mass measurement. For il- 
lustration we compare to distributions of two previously 
proposed variables. The first comparison is against the 
stransverse mass mT2 [3, HI , which is defined bj|f| 
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In Fig. 1(d) we show a distribution the uiT2{x = 0) 
distribution after applying the same basic dijet cuts used 
for Fig. 



1(a) 



The end-point of mn (x = GeV) distri- 
bution is also dependent on the physical quantity \pa\, 
and under the condition Q, 



™ m 2 X (x = 0) = 2 |po I 



(10) 



so the comparison is suitable. From the illustrative ar- 
rows it is clear that the end-point of the tot2 distribution 



is much more sensitive to the parent particle mass and 
that this end-point is in a region where Standard Model 
backgrounds arc suppressed. 



In Fig. 1(e) we show the distribution of another well- 



established SUSY mass-scale variable: the effective mass 
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where in this paper the sum is over the two jets with 
the largest transverse momenta. Since heavy particles 
are produced near threshold, the effective mass 14[ is 
expected to have an endpoint around 
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It can be seen from Fig. 1(c) that in contrast to mcT2, 
the m e ff endpoint is at a value for which the Standard 
Model background is small and that its characteristic 
value changes rapidly as a function of |j>o|- 

Even after the S cut designed to improve the peaked 
nature of the mcT2 endpoint the tot2 distribution shows 



more promise for the signal point investigated (Fig. 1(f) ). 



In conclusion, while it is true that the sharply peaked 
nature of the signal in the mcT2 distribution may lead to 
lower systematic errors in the fitting procedure, there are 
large backgrounds from Standard Model processes which 
also peak near this region. These backgrounds together 
with the weak dependence of the peak position on the 
physical parameter |po| suggest that mass measurement 
using this variable will be difficult. We suggest that for 
the decay process studied here, other existing mass mea- 
surement variables show more promising characteristics 
for SUSY mass measurement. 
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(a) mcT2(x = 300 GeV) distribution with just the dijet 
cuts applied. 
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(h) m^^ofv = 300 OeV) distribution after the ruts of |l 2| l . 
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(c) mcT2 {x = 300 GeV) distribution after the cuts of 112 
and an additional cut requiring \S\ < 30 GeV. 
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(d) mx2(x = 0) distribution with just the dijet cuts ap- 
plied. 
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(f) m.T2(x = 0) distribution after the cuts of 1 1 2D and an 
additional cut requiring <5| < 30 GeV. 



Figure 1: Results of the simulation described in the text for the SPSla SUSY benchmark point and 
various backgrounds. For this signal point the squark masses are in the range 500 < rriq < 600 GeV, 
the gluino mass is close to 600 GeV and the lightest neutralino mass is 97 GeV. The green (solid), blue 
(dotted) and red (dashed) arrows along the bottom of the plots show positions of the peaks ([5]) or 
end-points (fTOj) (fT2)) the distributions would be expected to have (under the condition ^) for processes 
with masses (my = 500 GeV, m x = 97 GeV), (m' Y = l.lmy, m x ) and (my = 2my, m x ) respectively. 
(The green and blue arrows lie almost on top of one another in the mcT2 plots). All plots correspond 
to integrated luminosity of lfb -1 . 
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